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CORRELATION OF AGN JET POWER WITH THE ENTROPY PROFILE 

IN COOLING FLOW CLUSTERS 

Fabio Pizzolato 1 , Tavish Kelly 2 , and Noam Soker 3 
ABSTRACT 

We find that the power of jets that inflate bubble pairs in cooling flow clusters of 
galaxies correlates with the size r a of the inner region where the entropy profile is flat, 
as well as with the gas mass in that region and the entropy floor— the entropy value at 
the center of the cluster. These correlations strengthen the cold feedback mechanism 
that is thought to operate in cooling flow clusters and during galaxy formation. In the 
cold feedback mechanism the central super-massive black hole (SMBH) is fed with cold 
clumps that originate in an extended region of the cooling flow volume, in particular 
from the inner region that has a flat entropy profile. Such a process ensures a tight 
feedback between radiative cooling and heating by the SMBH (the AGN). For a SMBH 
accretion efficiency (of converting mass to energy) of rjj = 0.1, we find the accretion 
rate to be M acc ~ 0.03(r a /10 kpc) L7 M yr _1 . This expression, as well as those for the 
gas mass and the entropy floor, although being crude, should be used instead of the 
Bondi accretion rate when studying AGN feedback. We find that the mass of molecular 
gas also correlates with the entropy profile parameters r a and M g , despite that the jet 
power does not correlate with the molecular gas mass. This further suggests that the 
entropy profile is a fundamental parameter determining cooling and feedback in cooling 
flow clusters. 



INTRODUCTION 



* 



It is now clear that in many clusters of galaxies moderate quantiti es of the intra-cluster med ium 
(ICM) are cooling to low t emperatures (T <C 10 5 K; see reviews by lPeterson &: Fabianl 120061 and 
McNamara fe Nulsenl 120071 ) . Moderate implies here that the mass cooling rate to low tempera- 
tures is much lower than the cooling rate expected without heating, but it is much larger than 
the accretion rate onto the supermassive black hole (SMBH ) at the center of th e cluster. The 



cooling mas s is either formin g stars (e.g., lO'Dea et al.l J2008I ; iRaffertv et al.ll2008l ). forming cold 



clouds (e.g-. lEdge et al.ll2010l ). accreted by the SMBH, or is expelled back to the ICM and heated 
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when it is shocked. In this moderate cooling flow model (jSoker et al.l l200ll ; ISoker k Davidl 120031 ; 
Sokerl 120041 ) . the SMBH is fed by cold clumps originating in an extended region (r ~ 5 — 30 kpc 



of the cool i ng flow, in a process termed the cold fe edback mechanism (jPizzolato k Sokerl 120051 ; 
Sokerl 120061 ; iPizzolatd 120071 ; iPizzolato k Sokerl |2010| ) . The cold feedback mechanism overcomes 



some severe problem s encountered by feedback models that are based on a ccreting gas directly 



from the hot phase (jSoker 



2006 



Soker et al. 



2009: IPizzolato fc Sokerl 12010 ) . such as the prob 



lems in the Bondi accretion (IMcNamara et al.l |201C]) . The cold feedback mechanism is compat 



ible with observations that show that heating cannot completely offset co o ling fe.g.. IWise et al 



2004 



McNamara et al 



Salome et al. 



2004 



Clarke et al 



to include it in feedback simulations (JGaspari et al. 



2004 



201 



W 



Hicks k Mushotzkvl J2005I : iBregman et al.l I2006J ; 



120081 ; IWilman et al.ll2009l; Peterson k Fabianll2006l ). and preliminary steps are taken 



The properties and behavior of the clumps that cool to low temperatures were studied by u s 
in previous papers (JPizzolato k Sokerl 120051 ; ISokerl 120061 ; IPizzolato! 120071 : IPizzolato k Sokerl |2010| ) . 
The distribution of cold clumps is compl icated, and there is no wa y for us to determine it. For that, 
in this paper we use our recent results IPizzolato k Sokerl (J2010I ) that are summarized in section 
[21 to derive in section [3] a semi -empirical phenom e nolog ical expressions for the accretion rate onto 
the SMBH. The main result of IPizzolato k Sokerl (|2010l ) and the previous papers is the sensitivity 
of cold clump evolution to the entropy profile. Readers interested only in the phenomenological 
formulae for accretion can skip section [2] and go directly to section [3j Our short summary is in 
section |U 



2. BASIC PROPERTIES OF COOLING CLUMPS 



The cold feedback mechanism assumes that the ICM is populated by a widespread assembly 
of cold clumps. We define the clump's density contrast 

6 = (p> - P )/p, (1) 

where p' is the mass density of the clump, and p is the mass density of the surrounding ICM. 



As shown in IPizzolato k Sokerl (|2005l . |2010| ) the fate of a clump critically depends on its 
initial overdensity. (i) Stable clumps. They start with a small initial overdensity, and they are 
eventually stabilised, (ii) Relatively unstable clumps. Moderately dense clumps, whose density 
contrast decreases as they fall in, but they manage to reach the center, (iii) Absolutely unstable 
clumps. These have 5 > 5c, and their density increase monotonically. 



The critical density is given by IPizzolato k Sokerl (J2010I ) 



[uj(8 c ,T)} 5 C = X, 



where 






T 



1 + 5 



(2) 
(3) 



depends on 5 and the ICM temperature T through the cooling function A. The left hand side of 
equation ([2]) is an implicit function of 5c, while the right hand side depends on the ICM properties 
and the radius of the clump according to 



X 



I C D 



g/a 



4 2 
W BV T cool 



where 



2 3 g d In K 

W BV 



(4) 



(5) 



5 r dlnr 

is the Brunt-Vaisala frequency, r coo i is the ICM cooling time, ais the clump's radius (assumed 
spherical) , g i s the gravitational acceleration, K is the entropy, and Cd — 0.75 is the drag coefficient 
(JKaiserl 12003). A plot of 5c as a function of x f° r several values of T is shown in Fig. [TJ Note 
that since x °t l/o, large clumps must be with a larger density contrast than small clumps to be 
unstable. 

From equation (JJ|) it is clear that the evolution of overdense, hence cooler, clumps is very 
sensitive to the cooling time and the entropy profile. Cooling of clumps is favored for a short 
cooling time and a flat entropy profile (by flat we refer to a moderate slope as well). The sensitivity 
of the clump evolution to the entropy profile is at heart of our present study. 



We define the entropy logarithmic gradient 



a 



dhiK 
dlnr 



(6) 



According to lDonahue et al.l (|2006l ) and lCavagnolo et al.l (|2009l ) in cooling flow clusters the entropy 
profile is flat, a < 1, in the center, and it has a ~ 1 in the outer regions. As in previous papers 
we take 10 < a < 100 pc. For 5c = 1 the condition for the clumps to be unstable is x ^ Xc — 0.1, 
where Xc is the critical value of x- Indeed, we define r a to be the radius where a = 0.5, and find 
that for our cluster sample 0.1 < x( r a) ^ 3. This condition on x translates to a condition on the 
cooling time, evaluated at ~ r a 

1 -1/2 



Tcool < 3 x 10 s 



10 kpcy VO.Olr 



-1/2 



{ 9 r) 1/2 
600 km s" 



a 

05 



Xc 1/2 yr- 



(7) 



It is interesting to note that iRaffertv et al.l (J2008I ) found that in clusters where the cooling time 
at a radius of 12 kpc is r coo i(12 kpc) < 8.5 x 10 8 yr = r s , high rate of star formation occurs. Equation 
([7]) above gives a limit of ~ 0.3r s , but for clumps of size a = 10 _3 r (a = 10 pc at r = 10 kpc) 
it gives a time scale of r coo i ~ r s . We conclude that a high cooling rate of over-dense clumps in 
the outer regions of the cooling flow, where a ~ 0.5 — 1, can occur, according to the cold feedback 
mechanism , when the cooling t ime is r coo i(10 kpc) < 10 9 yr, in agreement with observations of star 
formation (jRafferty et al.ll2008l ). 



In the inner regions where aCl, more clumps can cool. However, this region contains less 
mass than the outer regions, and it is more prone to AGN heating. Although there are much to 



be done, the cold feedback mechanism can gen erally account for th e limit of t s ~ 10 9 yr for a 
substantial star formation to occur, as found bv lRaffertv et al.l (|2008l ). 



It is observe d that major AGN activity occurs at a general time intervals of ~ 10 8 yr, (e.g., 
Wise et al.ll2007l for Hydra A). In the cold feedback mechanism the time scale is determined mainly 
by regions from where large quantity of gas migh t be cooling. This occurs where the entropy profile 
becomes flat (moving inward), r a ~ 3 — 30 kpc ([Donahue et al.l 120061 ; ICavagnolo et al.ll2009l ). The 
free fall time from these regions is Tff ~ 10 7 — 10 8 yr. The duty cycle might take somewhat longer. 
We attribute the time intervals between major AGN eruptions to the time it takes unstable clumps 
in large numbers to fall from these regions onto the SMBH. 



Sokerl (120081 ) considered both the duty cycle and the radiative cooling time r coo i(12 kpc) < 



8.5 x 10 8 yr = r s condition for high star formation rate in the frame of the cold feedback mechanism. 
Sokerl (|2008l ) proposed the criterion that the feedback cycle period must be longer than the radiative 
cooling time of dense blobs for large quantities of gas to cool to low temperature. It is possible that 
the two conditions are required for star formation: that cold blobs in large number are unstable 
and cool to low temperature, and that the AGN heating duty cycle is longer than the cooling time. 



SMBH MASS ACCRETION RATE 



As usually done, we take the power in the jets to be 

Ejets =rjjMa CC c 2 . 



(8) 



We take the energy of the jets to be equal to the energy in bubble pairs of clusters. The cavity power 
is calculated assuming 4pV of energy per cavity, and the buoyancy time scale th for its formation . 
The energy and power of bubbles are taken from the list c ompiled by [ McNam ara et al.l (|2010l ). 



with complementary material from iBirzan et al.l (|2004l ) and iRaffertv et al. 



(200 



a). 



The power of 

the jets is listed in the fifth column of Table 1. The clusters listed in Table 1 are those that we 
could find both entropy profile and jet (bubble; cavity) power in the literature. On the one hand 
this is the minimum energy, as some of the jets' energy goes to excite shocks and sound waves. On 
the other hand, we think the bubble filling is mainly non-relativistic, such that their energy should 
be (5/2PV). Over all, the usage of £j e ts = P cav = ApV/t^ for the power of cavities is adequate 
for the present purpose. However, for the reasons listed above w e think that the uncerta inties in 
the power of jets are much larger than the uncertainties listed bv lMcNamara et al.l (|2010l ). In the 
present Letter we do not include the uncertainties in the analysis, which is equal to assuming an 
equal (large) uncertainties to all points. 

The en tropy profiles, as well as some other cluster properties are taken from the catalog 
compiled by lCavagnolo et al.l (|2009l ). From the entropy profile we find the radius r a , listed in the 
second column of Table 1, defined to be the radius where a = 0.5 (equation |6|) . Inward to r a the 



entropy profile is flat. We also define the gas mass parameter 

4:7T 

M g = —p(r a )rl (9) 

listed in the third column of Table 1. We use this parameter rather than the total gas mass inside 
r a because in many cases the large bubbles near the center distort the density profile, and we are 
after a simple relation. We also list values of x{ r a) for a = O.Olr and gr a = (600 km s -1 ) 2 in the 
fourth column of Table 1. 

We could not find a simple correlation between the jet power and x, nor with the cooling time 
r coo i, and nor with M g /r r . nni ■ A short cooling time of r coo i < 10 9 yr at r ~ 10 kpc is a condition for 



high star formation rate (jRafferty et al.ll2008l ). However, once this condition is met it is not clear 



if the cooling time directly determines the mass accretion rate onto the SMBH. 

In the cold feedback mechanism, the entropy profile seems to be a more fundamental parameter 
for mass accretion rate onto the SMBH, once a short cooling time is established. It was already 
found that high star formation rate and high Ha lumino sity require that the entropy floor in the 



cluster be K$ = kTn 2 ' 3 < 30 keVcm 2 ( Voit et al.ll2008l ). We look for a correlation rather than 



a limit. In Fig. [2] we plot the jets (cavity pair) power versus the size of the flat entropy region 
r a (upper panel), and the gas mass parameter M g (middle panel). A correlation, albeit not very 
strong, is clearly seen in all panels of Fig. [2j We look for the simplest correlation in the logarithmic 
plot, namely, a linear one. The correlations we find are 

logP cav (erg s- 1 ) ~ 1.65 log r Q (kpc) + 42.6, (10) 

with Pearson R coefficient of R 2 = 0.67, and 

logP cav (erg s- 1 ) ~ 0.58 log M g (M ) + 38.5, (11) 

with R 2 = 0.56. Although these are not tight correlations, when the large scatter that is expected 
from the temporarily variations in the AGN power and its influence on the ICM is considered, these 
correlations have a merit. 

Two comments are in place here. First, M g strongly depends on r a , M g oc r 3 , hence the 
two correlations are not independent of each other. Second, because of the large scatter and 
small number of objects, we did not perform any deeper statistical analysis. We only point to 
the existence of a correlation, and the potential of using it to estimate an average, over time and 
many clusters, SMBH mass accretion rate. In a forthcoming more detailed study we will try to 
incorporate elliptical galaxies to the correltion. 

With the aid of equation ([8]) , we cast the desired phenomenological formulae for the accretion 
rate on the SMBH in the form 

^acc ^ 0.03 (-£—) ' {^-Y l M Q yT-\ (12) 



iokp C y vo.i 



Table 1: DATA ON CLUSTERS 



Cluster 



(kpc) 



M g (r Q ) 
(10 8 M G 



X 



Q) 



P 

1 cav 

'10 42 erg s" 



M mol 

;io 8 m ; 



K 
(keV cm 2 ) 



A478 


6.3 


30 


1.32 


100 


25 a 


7.8 


A1795 


15.1 


190 


1.11 


160 


55 6 


19 


Perseus 


20.8 


500 


2.78 


150 


350 c 


19.4 


A2199 


6.2 


12 


0.24 


270 


2.7 d 


13.3 


A2052 


7.4 


5.9 


0.11 


150 




9.5 


A4059 


3.0 


1.7 


0.19 


96 




7.1 


Cygnus A 


15.3 


150 


0.52 


3900 


15 d 


23.6 


A2597 


12.2 


110 


1.86 


67 


45 a 


10.6 


Hydra A 


11.6 


77 


0.83 


2000 


ll a 


13.3 


Centaurus 


1.2 


0.27 


0.89 


7.4 




2.2 


RBS797 


25.2 


1180 


2.77 


1200 




20.9 


2A 0335+096 


8.9 


50 


3.19 


24 




7.1 


A133 


14.6 


83 


0.71 


620 




17.3 


A262 


5.1 


3.4 


0.14 


9.7 


9.30 d 


10.6 


M87 


1.8 


0.55 


0.20 


6.0 


0.08 e 


3.5 


HCG62 


2.7 


0.99 


0.92 


3.90 




3.4 


MKW3S 


18.5 


130 


0.35 


410 


5.4 d 


23.9 



Notes: The meanings of the quantities are as follows. r a is the radius inward to which the entropy 



profile is fiat, i.e., a < 0.5. The entropy profiles and entropy floor Kq are from ICavagnolo et al 
(|2009l ); The gas mass parameter M g (r a ) is defined in equation ©; x is according to equation (JJJ; 



The power of the cavity pair, P rav , is calculated by d ividing the cavity pa i r ener gy, APV ', by the 
buoyant time, and taken from iMcNamara et al.l (120101) andlRaffertv et al.l (120061). The mole cular 



mass Mmni sources: fa) Edg e (2001); (b) 
Salome fc Combesl feooah ; (e) iTan et all pOO 



Salome h Combes 



2004 



Salome et all (J2006I ): (d) 



and 



M a 



M n 



0.04 ( Z 9 
\ 10 10 M P 



© 



("Mi 



&Y M B „-. 



(13) 



Considering the large uncertainties in the physical parameters and the large scatter in the graphs, 
there is no point to give the numerical values in the above expressions to a higher accuracy. 

The accretion rate does not increase as fast as the gas mass parameter M g , because as we 
move to larger distances from the center the cooling time increases, and the gas supply becomes 
less efficient. However, equations (|12p and (|13p support the cold feedback mechanism in suggesting 
that mass accreted onto the central SMBH is drained from an extended region, particularly from 
where the entropy profile is flat. 



McNamara et al.l (|2010i ) find no correlation between the molecular mass in clusters and the 



cavity power. We anal yzed the 10 clusters that are in our sample and have molecular mass in 



McNamara et al. 



Salome et al. 



(bgld) (taken from feded boOll : ISalome fc Combeslbood : ISalome fc Combes 



2004 



20061 ; iTan et al.l J20Q8). We confirmed their finding that the cavity power has no 
correlation with the molecular mass (we get a Pearson R coefficient of R 2 = 0.13). We do find a 
convincing correlation between the molecular mass M mo \ and the entropy profile quantities r a and 
M g , as presented in Fig. [3l We again look for a simple linear relation in the log-log plots, and find 



logM mol (M ) = 2.41ogr Q (kpc) + 6.8, 
with Pearson R coefficient of R 2 = 0.65, and 



logM mol (M ) = 0.89 log M g (M Q ) +0.50, 



(14) 



(15) 



with R 2 = 0.71. 



The above finding is very interesting. We find correlations of the entropy profile properties with 
both the cavity power and molecular mass. However, there is no correlation between the molecular 
mass and the cavity power. This shows two things. (1) The situation is not si mple. There is a 
large scatter probably because the AGN activity substantially varies with time (jNipoti fc Binney 
2005). (2) The entropy profile seems to be the basic quantity determining the cooling of gas, both 
to form cold reservoir and to feed the SMBH. 



We also correlate P cav and M mo \ with the entropy floor Kq (from lCavagnolo et al.ll2009i ). which 
is the value of the entropy in the inner flat region. The correlation with the molecular mass is poor 
(lower panel of Fig. [3]). As shown in the lower panels of Fig. [2j the jets' (cavity) power clearly 
increases with the entropy floor. The linear fitting gives 



logP cav (erg s- 1 ) = 2.31ogK (keV cm 2 ) +41.7, 



(16) 



with R 2 = 0.67. This behavior is opposite to what one would expect from a naive cooling interpreta- 
tion. A naive expectation is that for a lower entropy the cooling is more efficient. However, a higher 
entropy floor goes along with a more extended (larger r a ) region; we find an almost proportionality 



relation of -K"o(keV cm 2 ) ~ l.lr Q (kpc). We find that the cavity power (and possibly molecular mass) 
increases with the size of the flat entropy region. This is expected in the cold feedback mechanism 
(section [2|) , as long as the cooli ng time and the entr opy floor are below their respective thresholds 
of Trn ni(12 kpc) < 8.5 x 10 8 yr (JRaffertv et alJbood ) and K = kTu' 2 ^ < 30 keVcm 2 |Voit et al 
20081 ) . respectively. 



SUMMARY 



Our main finding is the correlation of the jet power that went to inflate the bubble -P C av> with 
the size of the flat entropy region r a (Fig. [2]). (By flat profile we refer to a moderate profile with 
a < 0.5. ) The positive correlation between P cav and properties of the flat entropy region in the 
inner regions of clusters is expected in the cold feedback mechanism. However, at this point the 
cold feedback mechanism does not give the form of the correlation. We have tried the simplest 
correlation in the log- log plot, and derived equations (fTUj) . (fl~2j) . The correlation with the gas mass 
parameter defined in equation ([9]) is given in equations (fTT|) . ([13]) . 



There is a slight possibility that the correlation arises from an underlying effect that is not 
related to the feedback process. For example, a larger flat entropy region will cause jets to deposit 
more of their energy in the inner region and inflate large bubbles. However, the positive correlation 
of r a with the molecular gas (Fig. [3]), despite that the molecular gas does not correlate with the jet 
power, bring us to reject this possibility. Further support to this rejection is the positive correlation 
with the value of the entropy floor (the flat part at the center) , as discussed in the last paragraph 
of section [3l We rather attribute a fundamental role to the flat entropy profile in the inner region 
in determining the cooling of the hot gas: The larger the flat entropy region is, the larger is the 
clumps' draining volume. 



The relations (|12p and (|13p derived here, and a similar one that can be derived from equation 
(|16p . can be used for an estimate of the average mass accretion rate onto the SMBH. They should 
replace t he Bondi accretion that does not fit accretion onto the SMBH in the center of clusters of 
galaxies (jSokerlhoOfil : ISoker et allhoOflUPizzolato k Sokenhoiol : iMcNamara et aklhoioh . 
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Fig. 1. — Plot of the critical overdensity 5c, implicitly defined by Equation [2] as a function of the 
parameter \ (Equation [2]) for four values of the ambient temperature: T = 0.5 keV (solid line), 
T = 1.0 keV (dotted line), T = 2.0 keV (dashed line) and T = 4.0 keV (dot-dashed line). 
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Fig. 2. — The cavity power (see Table 1) versus the size of the flat entropy region r a (upper panel), 
the gas mass parameter M g as given in equation ([9l middle panel), and the entropy floor Kq (lower 
panel) . 
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Fig. 3. — The molecular mass versus r a (upper panel) and versus M g (middle panel). For sources 
of data see Table 1. 



